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Abstract 

We investigate the possibility to constrain the pion distribution amplitude from the 
7*7* — ► 7r transition. For a surprisingly large range in the two photon virtualities we find 
that the transition form factor is essentially independent of the distribution amplitude. 
This in turn entails a parameter-free prediction of QCD. The 7*7* — > n, r( form factors 
are also briefly discussed. We estimate that experimental studies might be feasible at the 
existing e + e~ experiments BaBar, Belle, and CLEO. 
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1 Introduction 



One of the simplest exclusive observables is the form factor F P ^*) for transitions from a real 
or virtual photon to a pseudoscalar meson P. Its behavior at large momentum transfer is 
determined by the expansion of a product of two electromagnetic currents about light-like 
distances [|l|. The form factor then factorizes into a hard scattering amplitude, which is 
known including the first-order perturbative QCD corrections [|, and a soft matrix element, 
parameterized by a process independent meson distribution amplitude 

For space-like momentum transfer the form factors F P (*) can be accessed in e + e~ — > e + e~P. 
The measurement by CLEO |4j] for quasi-real photons and P = it, i], 7/ has renewed the interest 
in these quantities, and many papers have been devoted to their theoretical analysis, e.g. ||- 
HJ, to name a few. The CLEO data are consistent with distribution amplitudes of the pion, 



the 77, and the rf which are rather close to the asymptotic form, 

*As(o = |(i-e a ), a) 

where £ = 2x — 1, and x is the usual momentum fraction carried by the quark inside the meson. 

The purpose of the present article is to investigate the information contained in the form 
factors for 7*7* — > P transitions, beyond what we have already learned from the real-photon 
case. In Sect. ^| we cast the leading twist, next-to- leading order result for the 7*7* — > n form 
factor in a form useful for the purpose of our study. In Sect. ^|we attempt a critical appraisal 
of what we can and what we cannot deduce from the existing data on the transition 7*7 — > n. 
The following two sections explore the 7*- tt transition form factor in two different kinematical 
regimes. In Sect. ^ we briefly point out the specifics of the transitions to rj and rf mesons. 
Estimates of cross sections at the running experiments BaBar, Belle, and CLEO are given in 
Sect. 0, and we conclude in Sect. §. Some technical details concerning the a s corrections to 
F ny * are given in an appendix. 



2 The 7*— 7T transition form factor to leading twist 



Let us begin with the discussion of the 7*7* — > 7r form factor. The 7*7*7r vertex is parameter- 
ized by 

IV = -ie 2 F^j* (Q, uj) q a q' p , (2) 

where we use the convention eoi23 = 1- Here q and q' respectively denote the photon momenta 
corresponding to the Lorentz indices \i and v. We introduce the spacelike photon virtualities 
Q 2 = —q 2 , Q' 2 = —q' 2 , as well as 
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(Q 2 + Q 



'2\ 



Q 2 -Q 



(3) 



2^ ' 11 ~ Q 2 + Q' 2 ' 

The values of u range from —1 to 1, but due to Bose symmetry the transition form factor is 
symmetric in this variable: F T7 *(Q,u) = F irf »(Q,—ijj). 

To leading- twist accuracy, i.e., in the collinear approximation and using only the valence 
Fock state of the pion, the transition form factor F ny * reads |fj, |3j 



F nr (Q,u) 



3V2Q 
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(4) 
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Figure 1: Lowest order Feynman graph for the 7*7^) 
obtained by interchanging the photon vertices. 



7r transition. A second graph is 



The Feynman graphs contributing to leading order (LO) are shown in Fig. [TJ, and the next-to 
leading order (NLO) kernel JC(u,£,Q / ' fip) in the MS scheme is given in the appendix, fip 
and fin respectively denote the factorization and renormalization scales, both to be taken of 
order Q, and f n ~ 131 MeV is the well-known pion decay constant. <& n is the pion distribution 
amplitude, which we expand upon Gegenbauer polynomials C%/ 2 (£), the eigenfunctions of the 
leading-order evolution kernel for mesons ffl, 



E B n {f, F )Cl'^) 



n=2,4,. 



(5) 



The Gegenbauer coefficients are scale dependent, to LO they evolve according to 



B n (fii 



a. 



7n//?0 



(6) 



where fio is the starting scale of evolution and /3q = 11 — 2n//3. The anomalous dimensions 
7„ are positive numbers increasing with n, so that for hifip — > 00 any distribution amplitude 
evolves into $as with higher order terms in (gj) gradually becoming suppressed. Hence, the 
limiting behavior of the transition form factor for real photons is 



F nl (Q,u = ±l) 



V2Q : 



(7) 



which is a parameter-free QCD prediction once f n is known. Note, however, that this limit 
is only approached logarithmically, and that the anomalous dimensions are not very large at 
small n. For rif = 4 flavors one has 72/A) = 2/3, 74/A) ~ 0.97, •Jq/Po ~ 1.17. As n becomes 
large the 7„ grow logarithmically, 7 n ~ y hi(n + 1) being a good approximation already for 
n = 2. 

To NLO accuracy, the C^ 2 (£) are no longer eigenfunctions of the evolution, so that their 
coefficients do no evolve independently. Namely, B n (fip) at fip > fio depends on all coeffi- 
cients B 2 (fio), B n (fio). NLO evolution resums logarithms a 2 s log(fiF / fio) , compared with 
a s log(fiF / fJ>o) m LO evolution, and its effects will be more important when one evolves over a 
large interval in fip or when a s at the starting scale fio is large. 
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Using the Gegenbauer expansion (|), the integral in (f|) can be worked out analytically 
order by order in n, provided that lir is chosen to be independent of £. This results in 



n=2,4,. 



"n(u,/J, R ,Q/ 'l! F ) B n (fi F ] 



(8) 



where the lowest-order coefficients c n read 
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The analytical expressions of the a s corrections /C 2 and /C4 are rather lengthy and we refrain 
from showing them explicitly. Their dependence on ln(Q/fi F ) is partially compensated in 
F n ^*(Q,u) by the \i F dependence of the B n (n F ). Notice that no such compensation takes 
place for the [ar dependence to the order in a s we are working in. Unless stated otherwise we 
will in the following take li f = lir = Q, which is the virtuality of the quark propagators in 
Fig. [I| at £ = 0. The coefficients c n then depend weakly on Q via a s (Q). Other scale choices 
lead, as usual, to results differing by terms of 0(a 2 ), which is beyond the accuracy of our 
analysis. We finally remark that for fixed u the F nl * form factor only falls off like Q~ 2 at large 
Q, in contrast to the Q~ 2 Q'~ 2 oc Q~ A behavior of the vector meson dominance model fl~4|| . 

In order to visualize the sensitivity of the form factor to the Gegenbauer coefficients we plot 
the lowest coefficients c n (u) in Fig. [|. Here and in the following we use the two-loop expression 
of a s for n f = 4 flavors and = 305 MeV pf . We see a surprising behavior of the coefficients 
in uj. In the real-photon limit u — > 1 the form factor is sensitive to all Gegenbauer moments 
with approximately equal weight, c n (cu — 1) « 1. As soon however as one departs from this 
limit, the coefficients c n (uS) decrease and become ordered as c > c 2 > c 4 > . . . . Except for 
c this decrease is rather fast. If uj < 0.8, for instance, the second coefficient c 2 is less than 
40% of Co- For a wide range of uj, the form factor is essentially independent of the B n , unless 
they are unexpectedly large. The discussion of what can be learned from the measurement of 
F^rf naturally falls into two parts, concerning the kinematic regions Q' 2 -C Q 2 and Q' 2 ~ Q 2 , 
respectively. 



3 The real-photon limit 

Before embarking on the study of F 7r7 * let us discuss some issues in the analysis of the real- 
photon limit. At uj — 1 one has 



c n (u = 1) = 1 + as ^ R K, n (u = 1, Qfn F ) , 
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(10) 
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Figure 2: The coefficients c n (u) in the expansion (H) of the 7*-7r transition form factor. NLO 
corrections are included with [if = = Q, which is taken as 2 GeV. 



where the first few coefficients explicitly read 

3' 3 \72 6 n\ j ' 3 V 4500 75 ^ J K J 

The 7- 7r transition form factor thus approximately probes the sum 1 + J2n B n of Gegenbauer 
coefficients. Due to evolution and the running of a s the relative weights of the _B n (/io) in 
vary with Q, but only logarithmically. Extracting information on the Gegenbauer coefficients 
beyond their sum hence requires analyzing the form factor in a sufficiently large range of Q, say 
Qmin < Q < Qmax- While Q m ax is usually set by the available data, the choice of Q m in is an 
important source of theoretical uncertainty, since at lower Q power corrections to the leading 
twist result can become increasingly important. There are various approaches to describe such 
power corrections, all of which involve further assumptions or parameters. 

A different source of uncertainty is that in practice analyses of F nl data are performed with 
a truncation of the Gegenbauer series, setting all B n = for n> uq with some n . This may 
be seen as an analog to the determination of parton distributions, where one typically chooses 
a functional form of the parton distributions at the starting scale of evolution and then fits its 
parameters to data on inclusive processes. Notice that setting B n to zero for all values of //p is 
strictly speaking not consistent with NLO evolution, which generates nonzero B n {np) even if 
Bn(fJ'o) — 0. Miiller |TB[ has shown that this can lead to important effects on the shape of the 
distribution amplitude, especially in the endpoint regions and when a s at the starting scale 
is large. In the analysis M of the CLEO data, the impact of NLO evolution on the quantity 
F wl was however found to be small compared with the NLO corrections to the hard scattering 
kernel. In the numerical studies of the present work we will be concerned with values of fip = Q 
between 1 and 2 GeV. We choose /x = 1 GeV as the starting scale of evolution and take the 
LO formula @, using however the two-loop expression of a s . Since the purpose of this paper 
is not a precise determination of <&„-, this should be sufficiently accurate. At the same time it 
keeps the analysis procedure simple as it allows us to work with a finite number of Gegenbauer 
coefficients at all scales. 

The simplest analyses of the F 7r7 data assume no = 4. This leads to stable results on B2, 
although its value is subject to the uncertainties just discussed. Thus, for instance, Ref. 
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obtained the value B 2 (fi = 1 GeV) = —0.15 ± 0.04 to NLO accuracy in the MS scheme, using 
Hf = = Q and Q 2 min = 3 GeV 2 . Changing to fxp = Hr = Q as we prefer here, one finds 
5 2 (>o = lGeV) = -0.06 ± 0.03 for both Q 2 min = 2 GeV 2 and 3 GeV 2 within errors. This 
shift in the value of B 2 may be taken as an indication of the uncertainties due to uncalculated 
higher orders in the perturbative expansion. Melic et al. [13] have calculated the part of the 
0(a 2 ) corrections that allows one to determine the BLM scale for F wry . Taking the distribution 
amplitude $as they find [1r ~ Q/3 in the MS scheme. The corresponding NLO prediction for 
F nl is slightly below the CLEO data, in contrast to the one for <3>as an d Hr = Q, which is 



slightly above. Brodsky et al. |TTJ find consistency with the CLEO data when taking $as an d 
a yet lower renormalization scale together with a prescription for the saturation of a s , thus 
including effects beyond a leading twist perturbative analysis. 

If one allows for B 2 and B 4 in the analysis there is no unique result for the individual 
coefficients. Rather there is a linear correlation between B 2 and B 4 , which using (|6]), (Q), ( |TT]) 
one can approximate as 



(74-72)/% 

HQ) = /?,(//„) (i ^^(^4-^)1 l^g^l B 4 (fi ) + O(a 2 s 

\ 0.3 



^ Vc* - ic 2 )) ( ^4) 

« S 2 (/i ) + (l + 0.66a s (g)) f^y) B,(fi ) + 0(a 2 s ), (12) 

where we have taken fip = fin = Q and rif = 4 when going from the first to the second line. 
The function f(Q) includes the data on F nl and the term with cq in (|8|). It may also absorb 
possible power corrections as far as they are independent of Only the mild logarithmic Q 
dependence due to evolution and the running of a s restricts the allowed values of B 2 and -E> 4 to 
a finite region in parameter space. The factor multiplying 5 4 (/i ) in ( |T2"D varies between 1.31 
and 1.05 for Q 2 between 1 and 4 GeV 2 , the latter corresponding to the center of the highest 
Q 2 bin in the CLEO measurement. 

The experimental errors of the form factor data allow deviations from the linear correla- 
tion (|i~2D. This correlation is nicely illustrated by the x 2 - con t° ur plots in Fig. generated 
by MINUIT for our leading-twist NLO analysis. Comparison of the plots reveals that for 
Qmin = 3 GeV 2 the allowed parameter regions are enlarged since both the number of data 
points included in the fit and the lever arm in Q are smaller. On the other hand, power 
corrections should be less important in this case. We also observe from the figure that the 
la range for B 2 and £> 4 obtained from the fit with Q 2 min = 2 GeV 2 is embedded in that for 
Qmin = 3 GeV 2 , indicating that the fits are consistent with each other. Within the experi- 
mental errors, logarithmic effects suffice to describe the residual Q 2 dependence of the CLEO 
data for Q 2 F^iQ 2 ) above 2 GeV 2 . We emphasize that this finding does not prove that power 
corrections are indeed small in that region, it rather illustrates the difficulty to distinguish a 
power from a logarithmic behavior in Q 2 with data in the range between 2 and 8 GeV 2 . 

From the above exercise we conclude that only extreme values of |-B 2 | and I-B4I, say above 
1 or 2, are ruled out if higher order Gegenbauer coefficients are neglected. A compact way 
of presenting the result of the fit with two free coefficients is to use the linear combinations 
B 2 + -B4 and B 2 — B$, which have approximately uncorrelated errors. Taking Q 2 min = 2 GeV 2 
we obtain B 2 + B 4 = -0.06 ± 0.08 and B 2 - B 4 = 0.0 ± 0.9 at /x = 1 GeV. This illustrates 
that, within a leading twist NLO analysis, the CLEO data on the 7*7 — > it form factor is 
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Figure 3: 1 a and 2 a x 2 -contour plots for a two parameter fit to the CLEO data on F W1 
choosing Q 2 min = 2 GeV 2 (left) and 3 GeV 2 (right). Values of the Gegenbauer coefficients B 2 
and -B4 refer to a factorization scale of /io — 1 GeV. The data is fitted to the leading twist 
NLO expression @ with jj F = fx R = Q/y2. The y 2 is calculated from the errors on the data 
and does not include a theory error. 



insufficient for an unambiguous determination of the pion distribution amplitude, rather it 
approximately fixes the sum Yn B n to be close to zero. We cannot decide from the existing 
data on F ni whether the small value of the sum J2 n B n results from the cancellation of rather 
large individual terms or from the smallness of the B n themselves. Additional information 
from processes where $ ff enters in a different way, would be highly welcome to settle this issue. 
One candidate are the decays of charmonium states into pion pairs, where however the quality 
of the data as well as unsuppressed color octet contributions [|T7| do not permit a conclusive 
analysis at present. Many theoretical studies have been devoted to the elastic pion form factor, 



cf. e.g. []lT, lj| [3j|, ^TJ, but data at large Q 2 is scarce. A different way to probe the pion 
wave function is provided by diffractive dissociation of a pion into jets p2[ . It is however not 
clear at present to which precision information on Q n can be obtained in that process 



We will investigate below what data on the 7*- it transition form factor could contribute, 
but before this we should discuss the question of theoretical uncertainties in F n ^. The y 2 - 
contours in Fig. |3] only reflect the errors on the form factor data. We have already mentioned 
the uncertainty from the choice of factorization and renormalization scales. Much harder 
to estimate is the role of power corrections. Contributions from various dynamical sources 
have been studied in the literature, such as higher twist distribution amplitudes |24], ^5)1 , the 
nonperturbative behavior of a s (fi) in the infrared region HIT , 19| , or transverse momentum 
effects in the hard scattering subprocess |j[ ^ [LTj]. 



At this point we wish to comment on the analysis of the CLEO data presented in |f26 
Power corrections due to the hadronic component of the photon have been modeled there 
within QCD sum rules. The relative weights of B 2 (fio) and B^iiq) in F 7r7 then display a 
much stronger Q 2 dependence than in the leading twist case (|T2]), which leads to a much 
smaller allowed parameter region than in our Fig. [3|. A large part of the deviation between 
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the F %1 data and the result obtained with $as in that analysis is due to the inclusion of 
twist-four distribution amplitudes. Their shape is taken to be the asymptotic one, and for 



their normalization results from QCD sum rule calculations as given in [p4J] are used. No error 
is however assigned to this input in the analysis. Note that an uncertainty of about ±30% 
has been estimated for the relevant normalization constant 5 2 in |B7fl. Given this and the fact 



that only certain power corrections are taken into account in the QCD sum rule technique, 
we feel that the errors on the extracted Gegenbauer coefficients given in |26j are subject to a 
significant model dependence. 

Several of the analyses cited above find moderate but non-negligible power corrections to 
F W1 in the Q 2 range where most of the CLEO data is concentrated. We find that our limited 
ability to reliably describe, let alone to calculate power corrections precludes us from drawing 
firm quantitative conclusions and from discriminating between many of the available theory 
predictions for The best remedy to this situation we can see is data on at the highest 
possible Q 2 . Given the luminosity of the presently running experiments BaBar, Belle and 
CLEO, a substantial improvement should be possible over the statistics of the available data 
[[|, which is based on an integrated luminosity of about 3 fb _1 . 

Notice that the range of theory predictions on <&„. is considerable. For instance, the QCD 
sum rule analysis of Braun and Filyanov |28j gave -82(1 GeV) = 0.44 and -84(1 GeV) = 0.25 
under the assumption that higher order coefficients are negligible. On the other side, a still 
preliminary result from lattice QCD gives B2 = —0.41 ±0.06 at a low scale. Between these 
extremes many studies, using e.g. light-cone QCD sum rules with non-local condensates f30|, 



the transverse lattice [31], or the instanton model of QCD p2 obtain a distribution amplitude 
either slightly broader or slightly narrower than <3>as- 



4 The region Q' 2 < Q 2 

With the lessons from the real-photon limit in mind let us now investigate the region where 
one of the photons is slightly off-shell. In order to gain some insight in the importance of power 
corrections we estimate transverse momentum effects by employing the modified perturbative 
approach [[33], |2(J . It has been applied to the case of 7*- 7r transitions in || [7j , and to the case 
at hand in H. In this approach the expression (01) is replaced by 



F in *(Q,u) 



7T Z 



d 2 h -b, fi F ) K Q (y/l+fa Q b) exp [-S (f, b, Q, fi R 



(13) 



where Kq is the modified Bessel function of order zero, representing the Fourier transform 
of the leading-order hard scattering amplitude in momentum space. The quark-antiquark 
separation b is canonically conjugated to the usual transverse momentum kj_. The Sudakov 
exponent S describes gluonic radiative corrections not taken into account in the evolution of 
the wave function. For InQ — > 00 it suppresses all contributions to the integral except for 
those with small quark-antiquark separations. As a consequence the limiting behavior (0) 
emerges for 7- 7r transitions, as we have checked numerically by calculating expression (O) up 
to Q 2 = 10 15 GeV 2 . As b sets the interface between non-perturbative soft gluons contained in 
the hadronic wave function and perturbative soft gluon contributions resummed in the Sudakov 
factor, the factorization scale fip is taken as 1/6. For the renormalization scale we take the 
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Figure 4: Ratio of F nl * (Q, uj) in the modified perturbative approach and in the LO leading 
twist approximation at Q 2 = 4 GeV 2 (solid line), Q 2 = 2 GeV 2 (dashed line), and Q 2 = 1 GeV 2 
(dash-dotted line). Here we have used the wave function ( III] ) in the modified perturbative 
approach and the asymptotic pion distribution amplitude in the leading twist calculation. 



prescription hr = max {1/6, a/1 + 6^ Qi \A — <3} °f [13- Following [|4|, |2l] we take for 
the light-cone wave function in 6-space the simple form 



v/6 



$As(O ex P 



VT 2 / 2 
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in our estimate. Evaluation of the 7*7 — > 7r form factor in the modified perturbative approach 
using this wave function leads to very good agreement with the CLEO data @. 

In the kinematical range of interest here we find that the Sudakov factor only provides 
corrections of no more than 1.5% to F ni *, and it is thus good enough to retain only the 
leading logarithmic terms in S as given in In our kinematics the difference between the 
asymptotic result ([?]) and the expression (p~3|) is thus essentially due to the kj_-corrections to 
the hard scattering amplitude and not to the perturbative corrections contained in the Sudakov 
factor. In order to estimate the importance of power corrections we therefore compare fllTf ) 
with the leading twist result at LO rather than at NLO in a s . In Fig. [| we show the ratio 
between the form factor evaluated from ( |T4"D in the modified perturbative approach and the LO 
result calculated with $as m the leading-twist approximation. As we can see, the corrections 
are below 10% for values of Q 2 = 4 GeV 2 , but can go up to 30% for Q 2 = 1 GeV 2 . As expected, 
the importance of power corrections decreases as both photons become virtual. This is already 
signaled by the leading twist result. Indeed, the factor factor 1/(1 — i 2 uJ 2 ) in the convolution (f|) 
controls to which extent F ny * is sensitive to contributions from the end-point regions £ —>■ ±1, 
where the quark or antiquark in the pion becomes slow and soft effects can become important. 

Let us explore how much information on the pion distribution amplitude can be obtained 
from the uj dependence of F nl * under these circumstances. To this end, we plot in Fig. |5| the 
form factor for different choices of distribution amplitudes. These have been chosen to give the 
same value of J2 n B n — —0.06 at the scale fio = 1 GeV, so that up to small corrections they 
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Figure 5: The scaled form factor Q 2 F 7TJ *(Q,u) calculated to NLO in the leading-twist approx- 
imation at Q 2 = 4GeV 2 , using sample distribution amplitudes. The values of the Gegenbauer 
coefficients are quoted at the scale /i = 1 GeV. The curves are evaluated from B 2 = —0.06 
(solid line), from B 2 = -0.45, B 4 = 0.39 (dashed line) and from B 2 = 0.54, B 4 = -0.40 and 
Bq = —0.20 (dash-dotted line). All higher order Gegenbauer coefficients are taken to be zero. 



all give the same value of F nl in the NLO leading-twist analysis. We see that for Gegenbauer 
coefficients whose order of magnitude is not implausibly large compared to the theory estimates 
we have quoted above, one can obtain visible differences in the form factor. In part of the u 
range they can attain 15% and are thus marginally above the level where we have estimated 
that at Q 2 = 4 GeV 2 power corrections can make a reliable extraction of the B n problematic. 
While one clearly has not enough discriminating power to pin down individual coefficients B n , 
one can gain valid information beyond what can be inferred from real-photon data. Notably, 
one can check whether the small value of the sum Y^n Bn, to which F ni is mainly sensitive, 
results from the cancellation of rather large individual terms or from the smallness of the B n 
themselves. This type of evidence would be rather complementary to the quantitatively more 
precise information we have argued to be accessible from high Q 2 data on F W7 . 



As we have seen in Fig. ^ the contribution from Gegenbauer coefficients B n to F nl * becomes 
small as one goes away from u — 1, with a faster rate of decrease for increasing index n. To 
understand this, we observe that the hard scattering kernel in can be Taylor expanded 
around u = 0. Using the Gegenbauer expansion (|5[) of $7r(£)> we fi n d the relevant integral to 



5 The region Q' 2 ~ Q 2 



be 
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= £ u m / de(i-e a )r+ £ ^a-ar^ /2 (o (15) 

m=0,2,... J 1 n=2,4,... 1 

m=n,n+2,... 

for the LO contribution to F 7r7 » . The condition m > n in the sum involving the _B n incorporates 
the orthogonality of Gegenbauer polynomials, as the £ integrals for m < n are zero. We thus 
have the remarkable result that in the limit w-*0a Gegenbauer coefficient B n is suppressed 
in F %1 * by a power o; n . This also holds to NLO, as we will show in the appendix, and explains 
the behavior of the coefficients c n in Fig. |2|. Keeping only terms up to 0(u 2 ), we have 



F nr (Q,uj) = — =| 



1 _ 
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+ 0(u\a 2 ). (16) 



The limiting behavior for lu — > and to leading order in a s has been derived long ago in [35]. 
The a s -correction to the leading term, 



F nr (Q,uj) = y^JL 



1 _ ^1 

7T 



+ 0(u 2 ,a 2 s ) (17) 



has already been given in ||. 

Given the small numerical coefficients in front of u 2 , the u independent term in Eq. (|16|) 
dominates over a rather large range of u. Even at u ~ 0.6 the u 2 corrections amount to less 
than 15% if \B 2 \ < 0.5. Since higher coefficients B n are suppressed even more strongly, we 
conclude that in this range of uo the 7*-7r transition form factor is essentially flat in uj and 
independent of the pion distribution amplitude ^ 7T . To illustrate the quality of the small-a; 
approximations we compare in Fig. |6| the full result (Q) for F nl * with the expressions ([16]) 
and (|17D at Q = 2 GeV for a sample distribution amplitude given by B 2 = 0.54, B4 = —0.40, 
Bq = —0.20 at yU — 1 GeV. The full calculation is in agreement with the CLEO data for uj — > 1. 
We see that, although B 2 in our example is quite large and positive, both approximations are 
indeed very good for u < 0.6. 

We thus have a parameter-free prediction of QCD to leading-twist accuracy, which well 
deserves experimental verification. Any observed deviation from ([H]) beyond what can reason- 
ably be ascribed to 0(a 2 ) terms would be an unambiguous signal for power corrections. Only 
if the lowest Gegenbauer coefficients B n were extremely large would this conclusion become 
invalidated, but as we discussed in the previous section, such a scenario could already be ruled 
out using the region 10 ~ 1. We remark that according to our estimate in Sect. power 
corrections need not be negligibly small at moderate values of Q, even for u = 0. 

For small u, the relation (|17D has a status comparable to the famous expression of the cross 
section ratio R = a(e + e~ — > hadrons)/o-(e + e~ — ► fi + {i~). An important difference is however 
that the hard scale in R is timelike and requires one to stay out of the resonance mass region. 
In contrast, F W1 » involves spacelike virtualities and thus offers the possibility to explore the 
quality of a leading twist approximation down to moderate values of the hard scale. In that 
respect it is similar to a number of sum rules in inclusive deep inelastic scattering. The fact 
that (0) should hold in a wide region of to raises hope for experimental feasibility of this 
prediction. 
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Figure 6: Comparison of a NLO leading twist calculation of the scaled form factor Q 2 F nl * (Q, u) 
(solid line) with the approximations (|16D (dashed line) and ( |17D (dash-dotted line). The form 
factor is evaluated with fip — Hr — Q at Q = 2 GeV for a sample distribution amplitude with 
B 2 = 0.54, £? 4 = —0.40, Bq = —0.20, and B n = for n > 8. Note the suppressed zero on the 
y axis. 



One may ponder on whether at large Q 2 the form factor F nl * has the potential for a 
determination of a s . Our arguments in the appendix suggest that the suppression of B n 
by u n holds to all orders in perturbation theory. Higher orders a s coefficients in (|T7p for 
$,r = $as have been obtained in using the conformal operator product expansion. Given 



phenomeno logical or theoretical input on B n , the u n terms could at least be estimated, as 
could be the size of power corrections. Experimentally, the measurement of F nj * should be 
quite clean. We will however see that truly large Q 2 are not attainable at present facilities 
because of rather small cross sections. 



6 The 7*— r] and 7*— 7/ transition form factors 



Let us now discuss the 7*-P transition form factors for P — rj, rj'. For the valence Fock states 
of the mesons we choose a basis with the SU(3)p singlet \qq^) = \uu + dd + ss)/v3 and 
octet \qq^) = \uu + dd — 2ss) / \/6 , and the two-gluon state \gg). This has the advantage that 
the corresponding distribution amplitudes and mix under evolution, but evolves 
independently. The solution of the LO evolution equation for the octet distribution amplitude 
is given by (|5|) with Gegenbauer coefficients Bp n , whereas for the quark singlet and gluon one 
can write |]36| 
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Figure 7: Sample NLO Feynman graphs for the 7*7^*) — > P transition. 



with 



?(-)/ 



Pn 







(+) r(+) 



(19) 



where the Bp^ evolve as in (|6]) with positive anomalous dimensions 7^; 



(±) 



that conflicting results on 7^ and are found in the literature [56 



We remark in passing 
Notice that the two- 

gluon distribution amplitude vanishes in the asymptotic limit In /ip — > 00. It contributes 
to Pp 7 * only to order a s through the box graph shown in Fig. [?[ The corresponding hard 
scattering amplitude can be obtained by crossing from the NLO corrections to the Compton 
amplitude 7*p — > 7*p, which can be found in |57[] . 

In full analogy to the case of the 7r we can write the 7*-P transition form factor as a super- 
position of Gegenbauer coefficients with to dependent weights. In the limit u — *■ 1 one probes 



essentially the sums J2 n Bp! an d 2n Ppn °f coefficients appearing in the quark distribution 

. Using the mea- 



>(i) 



amplitudes and <3>p > . The real-photon limit has been analysed in |5|, 
surements of CLEO [[§] and L3 
the data on F, 



it was found in Ref. |38[ that within experimental accuracy 



and P ?? / 7 is compatible with the asymptotic forms of the quark distribution 



amplitudes and correspondingly vanishing gluon ones. 

When u} moves away from 1 the form factors become increasingly less sensitive to the higher 
order coefficients. As in the case of quarks, we find that the Gegenbauer coefficient Bpl or the 
gluon distribution amplitude first appears at order u n , despite the fact that is expanded 



upon Gegenbauer polynomials instead of C^ /z . In analogy to (|I7D we then obtain the 



3/2 



prediction 



(Q,u) 



yj/g 
3Q 2 



1 _ ^1 

71 



(20) 



The effective (process-dependent) decay constants are 



reft 

Jp 



/p 8) + 2V2£ 



(21) 



where 



(0\J$\P<P)) = if!Pp l » 



1. 



(22) 



are matrix elements of the SU(3)f singlet or octet axial vector currents. 
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The decomposition of the mesons states \r)) and |r/) on \qq^), \qq^ 8 '), \gg) and higher 
Fock states is completely general. It does not refer to 77 — 7/ mixing, which is known to occur 
empirically, and which relates the respective Fock state coefficients for the two mesons. In 
[[|(], f41| a quark-flavor mixing scheme has been proposed, which successfully describes many 
physical processes. In this scheme the physical mesons are obtained from two basis states, \rj q ) 
and \r] s ) by a unitary transformation with a mixing angle ip. The quark valence Fock states of 
\r] q ) and \r] s ) respectively are \uu + dd) / '\J2 and \ss). This mixing scheme can be justified to 
the extent that the distribution amplitudes are close to the asymptotic form; evolution does 



then practically not spoil this scheme. The decay constants (22) are given here by the mixing 
angle and the two basic decay constants, f q and f s . For the effective decay constants above 
one then finds [HOI 



re 

J lj 



cff 



5f q cos <p - y/2f s sin ip 



,es 5/ g siny? + y/2 f s cos <p 

Jri> — o ■ \ Z6 I 



The phenomenological values of the mixing parameters derived in [40| numerically give 



poff 



0.98/^ and = 1.62 f^. Similar results are obtained with the large N c values of the fp ||42| . 
At small u and large enough Q 2 the ratio of the 7*- 77 and r f—r( form factors constitutes an 
accurate measure of the effective decay constants. This can be used for a severe test of the 77—77' 



mixing scheme. As in the pion case, substantial deviations from the small-cu predictions (20) 
only occur if there are unexpectedly large power corrections or extremely large Gegenbauer 
coefficients. 

One may extend our analysis to the gluonic transitions g*g* — > 77, 77' in the same fashion 
as we discussed the electromagnetic ones. The case of the 77' has recently been investigated 



by Ali and Parkhomenko j43|, and we can follow their analysis closely. For u — > 1 the form 
factors not only depend on the Gegenbauer coefficients B P ' n but are also sensitive to variations 

of the Bpn, since in contrast to the electromagnetic case both quark and gluon distribution 
amplitudes now contribute to LO in a s . On the other hand, we can still apply our arguments 
for the suppression of the Gegenbauer coefficients (and hence of the full contribution from Qp 1 ^) 
in the limit ui — > 0. We find for the transition form factors 



F Pg , (Q, ui) = -jqI C p + 0(u 2 , a 2 s ) (24) 

with effective decay constants 

C v = V2f q cos ip - f s sin ip ps 0.32/,, , 

C r / = >/2/ g sin^ + / a cos^wl.99/ ff . (25) 

The ratio of the two form factors at small u is given by C v / C v > = — tan 9\ which, if measurable, 
would give access to the badly constrained 77 — 77' mixing parameter 9\ |41| . This mixing angle 
determines the relative decay strength of the 77 and 77' through a weak SU(3)f singlet current. 
Present estimates of the angle 9\ range from —2° to —10° ||41]| . 
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E e + 


E e - 


<^min 


Qimax 


i9 ■ 

u mm 


^max 


BaBar 


3.1 


9.0 


478 


284 


667 


275 


Belle 


3.5 


8.0 


154 


112 


524 


294 


CLEO 


5.3 


5.3 


227 


227 


314 


314 



Table 1: Minimal and maximal values of polar angles in the laboratory frame imposed by our 
cuts, as explained in the text. We also give the beam energies for the different experiments. 
Angles refer to the positron beam axis and are given in mrad, energies are given in GeV. 

7 Cross section estimates 

In order to assess the possibilities of the running experiments BaBar, Belle, and CLEO to 
investigate F ny * we will now estimate cross sections and see how they are affected by acceptance 
cuts. We first remark that at large Q 2 the cross section for e + e~ — > e + e~ n° depends on the 
two photon virtualities roughly like 



dQ 2 dQ' 2 Q 2 Q' 2 (Q 2 + Q' 2 ) 2 ' 

where the factors Q 2 and Q' 2 are due to the photon flux, and (Q 2 + Q' 2 ) 2 comes from the 
behavior of -F 7r7 », which as we have seen behaves like l/(Q 2 + Q' 2 ) with only a mild dependence 
on u. 

For acceptance cuts we will consider two scenarios. One is that both scattered leptons and 
the 7T° are seen in the detector. In kinematics where this is not possible because either Q 2 or Q' 2 
is too small to ensure a sufficiently large lepton scattering angle, one may envisage the detection 
only of one lepton and the 7r°. If experimental resolution permits, the four- momentum of the 
undetected lepton can then be reconstructed using four-momentum conservation. 

We define polar angles with respect to the direction of the positron beam, and take the 
convention that Q 2 corresponds to the photon radiated from the e + . As typical cuts we impose 
for tagged leptons a minimum energy of 1 GeV and transverse momentum of 100 MeV, and a 
minimum transverse momentum of 200 MeV for the pion. We further demand that a tagged e + 
has a polar angle a e + > a min in the detector. Likewise, we require a e - < a max for a tagged e~, 
and $ m i n < < $ max for the polar angle i?„. of the pion. For the minimal and maximal angles 
we take the values given in Table [IJ To estimate F^*^ we use the leading-twist NLO expression 
(fD, taking for simplicity the asymptotic form $as of the pion distribution amplitude. 

It turns out that the cuts with the most serious impact are the angular restrictions on the 
scattered leptons, while the precise values of the other cuts have only a mild influence. We 
show in Fig. |] the differential e + e~ cross section in the region Q 2 ~ Q' 2 . For Belle and CLEO 
we have imposed detection cuts for both scattered leptons. One clearly sees how for larger 
where at fixed Q 2 one of the photon virtualities becomes small, our cuts do have visible effects. 
Imposing the same cuts for BaBar leaves essentially no cross section in the kinematical region 
we are considering. This is due to the limited forward and backward coverage of the BaBar 
detector: for photon virtualities large enough to bring the scattered lepton within detector 
acceptance, the cross section is already minute because of its strong decrease (|26|) with Q 2 and 
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Figure 8: The differential cross section for e + e — > e + e it as a function of ui for Q 2 = 
2 GeV 2 (left) and Q 2 = 4GeV 2 (right). Notice that the scaling behavior fl26|) translates into 
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Figure 9: The differential cross section near lu = 1 for Q 2 = 2 GeV 2 (left) and Q 2 = 4 GeV 2 
(right). Corresponding curves with cuts for Belle are between those for BaBar and for CLEO. 
In this region, the behavior (f26|) becomes da / (d\og[l — lu] dQ 2 ) ~ Q~ 6 . 
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Q' 2 . The same holds if we allow the e + to be undetected. Only in the case where the e~ is 
untagged do we obtain a signal, not shown in Fig. |[ 

In Fig. we show the differential cross section in the region of u close to 1, requiring only 
the scattered e + and the n° to be observed. If instead we require the e~ and 7r° to be detected 
in the region of u close to —1, we obtain somewhat smaller cross sections in the case of Belle. 
For BaBar we have little change at Q 2 = 4GeV 2 but almost no signal left at Q 2 = 2GeV 2 . 
For CLEO with its symmetric geometry, there is of course no difference between the two cases. 

We have also investigated the production of an rj or rj', imposing the same cuts as described 
for the pion case. The corresponding cross sections scale approximately like the squared tran- 
sition form factors Fp 7 «, i.e., the mass differences between the tt, rj, and 7/ have only little 
effect on kinematics and phase space in the region we are investigating. 

Concerning the values of cross section estimated here, we recall the benchmark luminosity 
of 30 fb _1 per year of the B factories. With the numbers in Fig. [| and |9] we conclude that for 
Q 2 around 2 GeV 2 studies should be possible, both when u is around zero and when \u\ ~ 1. 
As Q 2 goes up to 4 GeV 2 and more, event statistics will become increasingly problematic, 
so that, unfortunately, we do not expect precision measurements to feasible in that region 
with the current experimental setups. Such studies would greatly benefit from high-luminosity 
upgrades of the .B-factories. 

We remark that our rate estimates here are restricted to \uj\ < 1 — 10~ 3 and thus do not 
include the real-photon limit. There the cross section will be much higher and, as emphasized 
at the end of Sect. §, measurements with good statistics should be possible at higher values of 

Q 2 . 

8 Summary 

We have investigated the 7*-7r transition form factor to leading-twist accuracy including a s 
corrections. The chief purpose of our analysis is to assess what can be learned about the pion 
distribution amplitude from experimental data on F 7r7 * . 

Our main idea is to use the expansion of ^ w on Gegenbauer polynomials, and to write F n ^* 
as a double series in the Gegenbauer coefficients B n and in powers of the variable u, which 
describes the difference of the two photon virtualities. We find that, contrary to what one may 
expect, it is very difficult to obtain information on the B n beyond what can be inferred from the 
case where one of the photons is quasi-real, which essentially constrains the sum of Gegenbauer 
coefficients. Only for u values close to but less than 1 can one get more information. Effects 
on F nl * of the order of 10% to 15% can be obtained using coefficients whose magnitude is not 
implausibly large compared to theory estimates. Data in that range of uo can, for instance, 
allow a check whether the small value of the sum J2 n B n , extracted from the CLEO real-photon 
data, results from cancellations of rather large individual terms or from the smallness of the B n 
themselves. This type of information would be a valuable input into other phenomenological 
studies of $„.. 

One of our main findings is that a Gegenbauer coefficient B n in $„. contributes to F 7r7 » with 
a weight proportional to uj n . For a large range of u, in fact for uj < 0.6, we find that the form 
factor is independent of the Gegenbauer coefficients to a high degree of accuracy. Although 
this is bad news for a determination of $> n , it entails a parameter free QCD prediction of the 
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7*7* — > vr form factor. Any clear deviation from this result observed in experiment would be 
an unambiguous signal for power corrections, provided that the lowest Gegenbauer coefficients 
B n are not extremely large, a scenario which could be ruled out using the region u ~ 1. In a 
wide region of u around zero, data on F nl * would thus permit one to test the quality of leading 
twist QCD in the spacelike region, in a similar fashion as sum rules in deep inelastic scattering. 

Analogous results hold for the 7*7* — > r\ and 7*7* — > rj' form factors. Note that here 
the gluon distribution amplitudes contribute indirectly through mixing in the evolution and 
directly to order a s . Their contribution is further suppressed for small u. Data on F^* and 
Fjf-y* in that region could rather cleanly determine of a linear combination of the flavor singlet 
and octet decay constants of these mesons, and thus give valuable information on 77-77' mixing. 

Cross section estimates of the process e + e _ — > e + e~n° at the running experiments BaBar, 
Belle, and CLEO indicate that it should be possible, although challenging, to measure the 
transition form factors for virtual photons up to about Q 2 ^4GeV 2 . Limiting factors for 
these measurements are luminosity and the acceptance for lepton detection in the forward and 
backward regions of the detector. Measurements with better statistics and at higher Q 2 might 
be feasible at high-luminosity upgrades of the present I?-factories. 

Concerning the real-photon limit, we argue that although the present data on Fp 1 favor a 
small value for the sum Y^n B n , more precise statements can only be made at the price of theory 
assumptions on the nature and size of power corrections to the leading- twist result. High statis- 
tics measurements at large virtualities Q 2 would greatly alleviate this problem and should be 
feasible at BaBar, Belle, and CLEO. They should be worthwhile, since the pseudoscalar meson 
distribution amplitudes are fundamental quantities describing meson structure and providing 
benchmark tests for nonperturbative methods in QCD. They are also an input required for the 
calculation of several phenomenologically important processes in hard-scattering approaches. 
An example are exclusive nonleptonic B meson decays into pseudoscalars, where a good un- 
derstanding of the strong interaction dynamics would enhance the prospects of extracting 
information on CP violation. 
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Appendix 

The NLO hard scattering kernel in (f|), evaluated in the MS scheme, reads 



1 r 



(1 + £u) ln(l - £u) + 4(1 



oj) ln(l 



6 



+ (1 + £u) ln 2 (l - £u) - (1 - u) ln 2 (l 
+ \ In % [2(1 + iui) ln(l - - 2(1 - 



u) ln(l - u) + 3(1 + £u) 



uj) - 9(1 + £u) 
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+ 



1 



6u; 2 (l -e 



2(1 + fu;)(l + & - 2 w 2 ) ln(l - £w) - 2(1 + + u - 2 w 2 ) ln(l - w) 



;i + - a; 2 ) ln 2 (l - £w) + (1 + w)(l - cj 2 ) ln 2 (l - w) 



a; 



3^(i-a 



In 



5! 

/4 



1 + £u) ln(l - £cj) - (1 + u) ln(l - u) 



(27) 



We will now show that the relevant convolution in the NLO part of F^* can be written as 

oo 

i+ e B n cr(® 



1 1 _£2 



e 2 ^ 



n=2,4,- 



E ^ m / ^(i-e 2 )pm(0+ E ^(i-apm(0c' /2 (0, (28) 



m=0,2 



n=2,4 
m=n,n+2 



where the p TO (£) are polynomials in £ of order m. Due to the orthogonality of the C^' 2 (£) the 
sum involving the B n is here again restricted to m > n. To show (j28l) it is enough to establish 
that /C(u;,£) can be expanded in a double Taylor series 



m=0 1=0 



E 4,^ 2fc K) 

fej=0 



2/ 



(29) 



i.e., that K is analytic in the two variables uj 2 and u 2 ^ 2 at uj 2 = u 2 ^ 2 = 0. One readily sees 
that K, is even in uj and £, which is a consequence of Bose symmetry and charge conjugation 
invariance. Further inspection of (]27f) shows that we can write 



jc = f(u 2 1 u 2 e) + 



u 2 e 



(30) 



where / and g are analytic in their variables around uj 2 = uj 2 ^ 2 = 0. Hence g is also analytic 
in uj 2 and uj 2 — uj 2 ^ 2 . Finally one can see from (^7|) that g is zero for £ = 1. The apparent 
pole at cj 2 — uj 2 £ 2 = thus cancels in (|30|), and g/{uJ 2 — uj 2 ^ 2 ) also has the required analyticity 
properties. 

It is not surprising that K, is analytic in the pair of variables uj and Apart from 
the appropriate subtraction of collinear singularities, the NLO hard scattering kernel is the 
amplitude for the partonic subprocess 7*(g) + 7*(?') — > + in the collinear limit, i.e., 
at (A; + k') 2 = 0. From the relations q' 2 = —Q 2 (l — uj) and (</ — A;') 2 = — Q 2 (l — u;£) we see 
that for Q 2 > analyticity in u; and u;£ around u = u;£ = is equivalent to analyticity in 
the spacelike invariants q' 2 and (q' — k') 2 . This suggests that our result ( Jl5| ) will generalize 
to higher orders in a s , provided appropriate analyticity properties of the collinear subtraction 
terms in the hard scattering kernel. 
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